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0\ ■ ABSTRACT 

o^ : 

^ . 

I We present the results of a pencil-beam survey of the Kuiper Belt using the 

1—5 I Keck 10-m telescope. A single 0.01 square degree field is imaged 29 times for 

'nI" I a total integration time of 4.8 hr. Combining exposures in software allows the 

I detection of Kuiper Belt Objects (KBOs) having visual magnitude my ^ 27.9. 

<N ; Two new KBOs are discovered. One object having my = 25.5 lies at a probable 

(N| I heliocentric distance R ~ 33 AU. The second object at my = 27.2 is located at 

^ I R ^ 44 AU. Both KBOs have diameters of about 50 km, assuming comet-like 

in ; albedos of 4%. 

^ ' Data from all surveys are pooled to construct the luminosity function from 

^ I rrifi = 20 to 27. The cumulative number of objects per square degree, S(< m^j), 

■ is fitted to a power law of the form log^Q E = ainiR — 23.5), where the slope 
)> i' 

I ' a = 0.52 ± 0.02. Differences between slopes reported in the literature are 

?-j ■ due mainly to which survey data are incorporated in the fit, and not to the 

^ ■ method of analysis. The luminosity function is consistent with a power-law size 

• distribution for objects having diameters s = 50-500 km; dN oc s"'' ds, where the 

■ differential size index q = 3.6 ± 0.1. The distribution is such that the smallest 
5^ ■ objects possess most of the surface area, but the largest bodies contain the bulk 

of the mass. We estimate to order-of-magnitude that 0.2M^ and 1 x lO^'' comet 
progenitors lie between 30 and 50 AU. Though our inferred size index nearly 
matches that derived by Dohnanyi (1969), it is unknown whether catastrophic 
collisions are responsible for shaping the size distribution. Impact strengths 
may increase strongly with size from 50 to 500 km, whereas the derivation by 
Dohnanyi (1969) assumes impact strength to be independent of size. In the 
present-day Belt, coUisional lifetimes of KBOs having diameters 50-500 km 
exceed the age of the Solar System by at least 2 orders of magnitude, assuming 
bodies consist of solid, cohesive rock. Implications of the absence of detections 
of classical KBOs beyond 50 AU are discussed. 



^Alfred P. Sloan Research Fellow 
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Subject headings: Kuiper Belt - comets: general - solar system: formation 

1. INTRODUCTION 

Beyond the orbit of Neptune lies a disk of remnant planetesimals known as the Kuiper 
Belt. As outhned by Jewitt, Luu, & Trujillo (1998, hereafter JLT98), the ~100 known 
Kuiper Belt Objects (KBOs) divide into 3 dynamical classes. 

1. Classical KBOs reside in low eccentricity, low inclination orbits beyond 40 AU 
(JLT98). They are not associated with mean motion resonances with Neptune. 

2. Resonant KBOs, of which Pluto is the largest known member, have orbital periods 
commensurate with that of Neptune and are protected against close encounters 
with that planet (Malhotra 1996). They possess moderately high eccentricities 
and inclinations, possibly excited by Neptune during a transient period of orbital 
migration (JLT98; Malhotra 1995; Malhotra, Duncan, & Levison 1999, hereafter 
MDL99; and references therein). 

3. Scattered KBOs, of which 1996 TLee is one member (Luu et al. 1997), occupy large, 
highly eccentric and inclined orbits, the result of close encounters with Neptune 
(Duncan & Levison 1997). 

Of the 3 populations, the classical Kuiper Belt appears the most untouched dynamically. 
However, the record of primordial conditions preserved by the classical Belt, as observed 
today out to 50 AU, seems heavily weathered. Recent surveys estimate the mass of the 
observable Kuiper Belt within 50 AU to be a few x O.IM^ (e.g., Luu & Jewitt 1998; this 
paper). Hamid, Marsden, & Whipple (1968) used the trajectories of short-period comets to 
set an upper mass limit of 1.3Mq on a smooth ring within 50 AU. The mass of the nearby 
Belt is therefore ~10-100 x smaller than the ~15M0 extrapolated from the condensable 
material of the outer giant planets. Mass depletion since the time of formation is also 
suggested by the existence of bodies as large as Pluto. To build bodies of this size at 36 AU 
within 100 Myr (the estimated formation time of Neptune), the standard model of pairwise 
planetesimal accretion requires the primordial disk to have at least ~1-10 from 29 to 
41 AU, depending on the assumed sizes of seed planetesimals (Kenyon & Luu 1998; cf. 
Stern & Colwell 1997). 
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The cause of the presumed depletion is unclear, but Neptune is considered a prime 
suspect. Duncan. Levison, & Budd (1995) calculate that Neptune, when fixed in its present 
orbit, can gravitationally eject more than 90% of the Belt mass inside 39 AU over the age 
of the Solar System. About 50-90% of the mass between 39 and 50 AU may be depleted 
by gravitational perturbations alone. These simulations are sensitive to assumed initial 
eccentricities and inclinations of test particles. Collisions are also proposed to explain the 
missing mass, either by nudging objects into unstable orbits (Davis & Farinella 1997), or 
by grinding bodies down to dust to be transported by radiation pressure (Stern & Colwell 
1997). In both cases, however, only the smallest KBOs may be significantly depleted. 
CoUisionally relaxed populations place most of their mass in the largest bodies, but 
most of their geometric cross section in the smallest members. Collisions might therefore 
preferentially grind down the smallest objects, leaving the largest bodies undisruptcd and 
the total mass mostly intact. This expectation is borne out in computations by Davis & 
Farinella (1997). Important caveats for all collisional simulations of the Kuiper Belt include 
oversimplified prescriptions for the impact strengths of KBOs, reflecting our ignorance 
of their internal structure. Erosive velocities are thought to be gravitationally stirred by 
Neptune within 50 AU, but physically motivated estimates of the velocity dispersion have 
yet to be made in these simulations. 

However large or small, the destructive influence of Neptune on the Kuiper Belt may be 
limited to within the location of its outermost 2:1 resonance at 48 AU. This idea has led to 
speculation that the surface density of Belt material rises by ~2 orders of magnitude to its 
assumed primordial value somewhere beyond this radius (Stern 1996; Stern & Colwell 1997; 
MDL99). However, no classical KBO beyond 50 AU has yet been discovered. Assuming the 
shape of the KBO size distribution does not change with distance, JLT98 flnd by Monte 
Carlo simulation that their observations arc consistent with an edge to the classical Kuiper 
Belt at 50 AU — a "Kuiper Cliff". The flrst theoretical constraints on the classical Belt mass 
beyond 50 AU are provided by Ward & Hahn (1998). They flnd under certain conditions 
that a Belt containing 1.6M© from 48 to 75 AU damps Neptune's eccentricity to its current 
observed value of 0.009 by the action of apsidal density waves. They calculate that the 
addition of 10 times more mass in this region (masses comparable to those expected in 
the minimum-mass solar nebula) would reduce Neptune's eccentricity to less than lO"^*^ 
over the age of the Solar System. In these computations, the outer Belt is assumed to 
consist predominantly of small bodies (diameters <C 140 km) so that velocity dispersions 
are sufficiently low to sustain wave action. 

In the absence of any direct observations of the Kuiper Belt beyond 50 AU, we 
undertook a pilot survey utilizing the Keck 10-m telescope. A single 600s exposure on Keck 
can achieve a depth mv ~ 26, allowing objects 100 km in diameter with comet-like albedos 
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to be seen out to distances just beyond 50 AU. Combining exposures in software enables 
the detection of such bodies inside 70 AU. Our primary aim was to constrain the KBO 
luminosity function out to my ^28; in this goal, we succeeded. Our principal hope was to 
directly image the Kuiper Belt for the first time beyond 50 AU; this wish remains to be 
fulfilled at present. 

Observations are described in §0. Methods of data reduction and search strategies are 
set forth in §|^. Results, including actual detections and our construction of the luminosity 
function from = 20 to 28, are presented in Implications of our results on the size, 
mass, and distance distributions of KBOs are discussed in Our principal findings are 
summarized in 

2. OBSERVATIONS 

Data were taken on 31 August 1997 UT using the Keck II 10-m telescope atop Mauna 
Kea in Hawaii. The Low-Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) was 
mounted at Cassegrain focus and employed in direct imaging mode. The plate scale on 
LRIS's Textronix CCD was 0.215 "/pixel. The camera had a useable (vignetted) field of 
view of 5.67 x 7.34 square arcminutes (1582 x 2048 square pixels = 0.0115 square degrees). 
A standard V filter was used. The choice of V over R was motivated by lower sky brightness 
and greater solar flux at V. While some KBOs have higher reflectances at R, others also 
appear neutral (Tegler & Romanishin 1998). 

We searched for KBOs in a single, relatively star-free field at opposition [a = 22'^54™54*, 
S = -6°20'34" (J2000)].i Twenty-nine exposures, each 600s in duration, were recorded 
of this field. Data were read out from the CCD through two amplifiers operating 
simultaneously; this procedure halved the readout time to 60s at the cost of introducing 
small differences in the amount of noise between chip halves. Each frame was offset in 
position by ~5-100" relative to other frames; our dithering routine enabled the construction 
of high-fidelity flatfields ("skyflats") from the science data themselves (see § |3.1| ). Provided 
the Keck telescope functioned properly, our duty cycle efficiency was nearly 90%. Unforeseen 
crashes in the mirror alignment software limited our total effective integration time to 4.8 
hr over a 6.2 hr baseline. 

A Landolt field (Landolt 1992) provided photometric standards. The seeing ranged from 
0'.'65 to I'.'O full- width at half- maximum (FWHM), with the median seeing equal to 0'.'75. 



^The field happened to be located 43° away in ecliptic longitude from Neptune. 
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3. DATA REDUCTION AND SEARCH STRATEGY 

Kuiper Belt Object candidates are identified by their parallax motions (of order "/hr) 
against the fixed stars. We employed two search methods: a simple blinking of individual 
frames to visually scan for slow-moving objects, and a deep, recombinative blinking 
approach which blends the search algorithms of Gladman et al. (1998) and Cochran et al. 
(1995). 

Observations of candidates over a single night are insufficient to constrain orbital 
parameters and to prove membership in the Kuiper Belt. Candidates might instead 
be eccentric, near-Earth asteroids whose apparent motions mimic those of true KBOs. 
However, as discussed by Luu & Jewitt (1998), the possibihty of mistaken identity appears 
remote, since masquerading slow-moving objects have not appeared in their many surveys 
to date. We proceed on the assumption that our (small) field is likewise uncontaminated. 



3.1. Shallow Survey: Basic Blinking 

All image processing described in this paper was performed with the Interactive Data 
Language (IDL) software package. The 29 science frames were first corrected for CCD bias 
and pixel-to-pixel variations in gain (flat-fielded). For each science frame, a tailored flatfleld 
was constructed from the median of the other 28 dithered science frames.! The fact that 
the images were dithered ensured that each CCD pixel sampled the flat sky several times. 

Each flattened frame had its mean sky value subtracted and its flux normalized by 
scaling eleven bright, unsaturated stars distributed across the entire frame. Position offsets 
required to align the dithered images were obtained by minimizing stellar residuals of 
frames subtracted pairwise. 

Aligned images were blinked and visually scanned for slow-moving objects. Three 
images of comparable seeing, spaced about 1 hr apart, were blinked per session. Four 
triplets were blinked in all, including the flrst and last frames of the night. Results of this 
comparatively shallow survey are presented in § [4.1| . 



Seven additional frames from other observations during the same night were included in the median 
flatfield. 
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3.2. Deep Survey: Forward-Reverse Recombinative Blinking 

The basic idea underlying our deep survey is simple. Images are stacked and shifted 
on top of each other according to a hypothetical KBO proper motion. The shifted stack 
of images is co-added to form a recombination image. While stationary objects appear 
smeared in the recombination image, an object whose motion matches that assumed has 
its signal strengthened and appears as a single seeing disk. Thus, all collected photons are 
used to identify KBOs too faint to rise above the noise of an individual image. 

To reduce confusion and noise in the recombination image, it is desirable to remove 
non-KBO sources of emission from individual frames before co-adding. Towards this end, 
we subtracted from each individual frame the median of the other 28 (aligned) frames.^ 
Extended, stationary, low surface brightness emission (from resolved galaxies, for example) 
was thereby mostly removed from individual images. Some pixel positions did not have the 
full overlap of all 29 frames because of our dithering routine; these were purged to ensure 
uniform statistics. 

Cosmic rays and asteroid streaks remained in the median-subtracted frames. 
Substantial residuals from stationary point sources were also left behind, a consequence 
of frame-to-frame seeing variations. All three non-KBO sources of emission were largely 
eliminated by clipping high-valued pixels from the shifted stack of images. After 
experimenting with various schemes, we decided to clip the 5 highest values from each 
column of 29 pixels and average the remaining 24 values.! Columns not having the full 
overlap of all 29 frames due to the shifting process were purged altogether. Finally, to the 
clipped mean image we added a positive constant frame to restore the average background 
level to zero. The resulting (rectangular) array constituted our recombination image, which 
appeared satisfyingly clean aside from a few well-localized and easily recognizable residuals 
from bloomed stars. 

The proper motion vector of a KBO is described by its amplitude, /x, and its apparent 
inclination angle, 9, relative to the ecliptic as seen on the CCD. Following Gladman et 
al. (1998), we visually searched for KBO candidates by blinking, in any one session, 
4 recombination images corresponding to 4 successively higher amplitudes along one 
inclination. Objects characteristically came into focus and then smeared as their actual 



"^linage processing for the deep survey began with the sky-subtracted, flux-normalized images from the 
shallow survey. 

^Taking the median of all 29 frames as an example of an alternative scheme generated a still noisier 
background than averaging 24 frames {amedian-29 / (^average^2'i ^ \J ^ '^■^^)- 
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rates of motion were approached and passed. Recombination amplitudes ranged from 
/i = 1.1 to 6.3 " /hi in steps of A/i = 0.4 "/hr. Inclinations ranged from 9 = —5 to 5° in 
steps of A^^ = 5°.i These ranges cover proper motions (as seen at opposition) of KBOs 
moving on prograde, circular, heliocentric orbits with semi-major axes R of 20-120 AU and 
actual inclinations i of up to 30°. i 

Roughly 130 artificial KBOs were implanted at random locations and searched for 
simultaneously with true KBO candidates. Their magnitudes were spread uniformly 
between my = 26 and 29, and their orbital parameters were chosen randomly within the 
ranges cited above. Their presence in recombination images trained the eye to recognize 
bona fide KBOs, and their rate of recovery provided an estimate of true KBO detection 
efficiency as a function of magnitude. Differences between artificial objects' given and 
recovered properties {6mv, Sn, 66) provided estimates of systematic errors in the parameters 
of true candidates. Recombination spacings were just small enough to detect artificially 
implanted KBOs in at least 2 recombination frames. 

A list of KBO candidates was made containing objects (including artificial ones) which 
(1) focussed and de-focussed in the correct manner, (2) appeared in at least 2 adjacent 
recombination images, (3) did not appear as a single hot pixel in any 1 image, and (4) were 
not situated too close to the noisy environs of stellar /asteroidal residuals. Objects in this 
list had their magnitudes and proper motions subsequently refined on a grid of resolution 
A/i = 0.1 "/hr and A6 = 1?25. This process involved extracting square subframes 20 pixels 
wide surrounding each candidate and recombining them on the finer grid. Simply selecting 
the grid point {n, 9) for which counts inside a circular sampling aperture were maximized 
proved too simplistic a procedure. Often the maximum-count image simply pushed hot 
noise pixels into our sampling aperture. In practice, we selected the best recombination 
image based on visual appearance, a well-behaved fiux profile, and in the few cases where 
we could not decide, maximal counts. 

Without a second night to confirm the reality of our candidates, visual surveys of this 
kind are more prone to false detections. Even apart from human bias, noisy pixels may still 
conspire to masquerade as slow-moving objects. To estimate the number of false detections 
in our candidate list, we repeated our entire deep search on images recombined in the 
reverse direction. Reverse in this case actually means in the apparent prograde direction, 
since proper motions of KBOs at opposition are dominated by the Earth's parallax motion 



^Negative (positive) implies motion near a descending (ascending) node. 

small correction term due to the fact that our field was 4?5 away from opposition in ecliptic longitude 
was included in calculating these ranges. 
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and must appear retrograde. Artificial, apparently prograde objects were also randomly 
inserted in individual frames and searched for in images recombined in the reverse direction. 
In this reverse survey, which suffered the same kinds of errors as afflicted the forward 
survey, objects that fulfilled the four requirements listed above and that turned out not 
to be artificially inserted were deemed chance alignments of noise. We refer to these as 
"reverse survey noise objects." Statistical confidence in the detection of real KBOs demands 
that the number of candidates detected in the forward survey exceed the number detected 
in the reverse survey plus the uncertainty in the latter number. In practice, we blinked 
recombination frames without knowing whether they were recombined in the forward or 
reverse directions, thereby avoiding another potential source of human bias. 



4. RESULTS 

4.1. Shallow Survey Results 

One KBO was discovered by blinking individual images. The object, hereafter OBJl, 
appears at the ~5.5cr level in 21 out of 29 frames. In the other 8 frames, light from 
OBJl had fallen off the CCD chip as a consequence of our dithering routine. Figure [1| 
displays our newly discovered object in 3 consecutive exposures, and Table [l| summarizes its 
measured and inferred properties. Its motion over 6.2 hr is consistent with being uniform; a 
best-fit line through centroid positions yields proper motion parameters /i = 3.83 "/hr and 
9 = — 1?0. The corresponding heliocentric distance and inclination for an assumed circular 
orbit are R = 32.9 AU and i = 4?5. With a measured visual magnitude of my = 25.5 ± 0.3 
{la dispersion among 21 measurements), the object is 56 ± 6 km in diameter, assuming it 
has a comet-like visual albedo of 0.04 (Allen 1973). 

The area of sky covered by our shallow survey was As = 0.0102 square degrees, after 
correcting approximately for dithering losses (-11%) and area taken up by bright stars and 



galaxies (-0.5%). A discussion of the cumulative luminosity function is reserved for §|0. 



4.2. Deep Survey Results 

4-2.1. Artificial Object Recovery 

After refining estimates of candidates' magnitudes and proper motions, we culled 
artificial objects from the candidate list. Figure |^ displays our recovery rate rj of artificial 
objects as a function of their given my, for both forward and reverse surveys. In both 



.1, * X'J 

K??!|.UT 07:-59 i :^>V-^'F0y 46"x27''-. z .■'v^j 



■■■■ .■ -^if^.-;-:; ■ • 




■■■■ - --iSm . ■ ■ 



I. ^ 

i ■ 
d ■ > 



■ ■ ■ ■ ■ ■ 




! 



-10:42 



JHi .1-1 --J ■■- ' 

, . V 4 ,^ — J" 



, . : ■ - ...V- . ;^:V%:- 



- mom ■ ■ ■ ■ . ™ 



I • 1 



;vsyL. 1r3;QS^ /i^K-:--: 

.t^J^ . - . .. . 



: . L "^ 4 ? -= " ■ ■ ■ 




■7 I- ■ 



Fig. 1. — Individual exposures of OBJl, with time and angular scale indicated. The object 
appears in a total of 21 frames. 
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surveys, the rate of recovery was similar, falling from 100% near my = 27.3 to 0% at 
my = 28.4. The datasets were combined and fitted to the function 



(Gladman et al. 1998). The fit yields a detection efficiency which falls to 50% at my (50%) 
= 27.94, over a characteristic width W = 0.38 mag. Satisfyingly, my (50%) is only 0.04 mag 
brighter than the nominal 3a limit obtained by reducing the noise of an individual image 
by 

Figure |^ plots 5my, 6fi, and 69 — differences between artificial objects' given and 
recovered properties — versus their given magnitude. To clarify possible trends with 
increasing magnitude, we also plot averages and standard deviations within bins of width 
0.5 mag; these points are positioned at the centers of each bin. 

We note first that there is no significant bias in our estimation of parameters; averages 
Sniv, S/j,, and 69 are consistent with being zero. The scatter, however, is significant. We 
adopt the scatter in 6my as our estimate of the uncertainty in true candidates' magnitudes; 
the la dispersion increases from 0.22 mag near my = 27 to 0.33 mag near my = 28. The 
analogous la uncertainty in /j, ranges from 0.05 to 0.17 "/hr, and the la uncertainty in 
9 ranges from 2?3 to 3?2. These results for /x and 9 appear reasonable. A difference of 
A/i = 0.1 "/hr over a time At = 6.2 hr smears images by A/i x At = 0'.'6 — about half the 
value of the worst seeing during our observations. A difference of A9 = 2?5 at a fixed, typical 
amplitude of /x = 3 "/hr smears images over a comparable distance: /i x A9 x At = 0"8. 



A second KBO, hereafter 0BJ2, was discovered by blinking recombination frames. 
Figure | presents the best recombination image of 0BJ2, surrounded by images of the same 
object recombined at adjacent points on the (/i, 9) grid. Its smearing pattern is identical 
to those of artificially planted objects having similar motions. Properties of 0BJ2 are 
summarized in Table |1|. Its visual magnitude is my = 27.22 ± 0.22, and its proper motion 
parameters are /i = 2.92 ± 0.05 "/hr and ^ = 0° ± 2?3. Alternative recombination frames for 
0BJ2 were constructed by clipping the top 1 pixel out of each column of 29 pixels and then 
averaging the remaining 28 values. Exactly the same parameters for 0BJ2 were obtained. 
If we assume a comet-like visual albedo of pv = 0.04, these measurements are consistent 
with those of an object 46 ± 6 km in diameter, occupying a circular, uninclined orbit of 




(1) 



4.2.2. True Object Discovery and Upper Limits 
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Fig. 2. — Rate of recovery of artificially implanted objects versus their given magnitude. 
Detection efficiencies from forward and reverse surveys are averaged and fitted to equation 
(0), shown ClS db solid line. 
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Fig. 3. — Differences between given and recovered properties of artificial objects. Solid circles 
indicate objects recovered in the forward survey, and crosses denote objects recovered in the 
reverse survey. Error bars reflect ±1(7 dispersions in bins of width 0.5 mag. Since average 
differences are consistent with being zero, we conclude that our measurements of my, /i, and 
9 are not biassed. We adopt the dispersions to be our measurement uncertainties. 
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semi-major axis R = 43.9 ± 0.8 AU. 

Confidence in the reality of OB J2 is furtlier bolstered by Figure in which we compare 
cumulative numbers of objects detected in forward and reverse surveys. No false alarm 
went off in the reverse survey at the magnitude of 0BJ2; the object distinguishes itself as 
the brightest detection at 5.5cr. 

By contrast, we view all candidates in the forward survey fainter than my = 27 A as 
false detections, partly because their numbers do not exceed those in the reverse survey. No 
object in both surveys is as visually convincing as 0BJ2; many other candidates vanished 
at several (but not all) adjacent recombination gridpoints. Moreover, regarding the last of 



the 4 search criteria set forth in § |3.2| , it was occasionally unclear when an object was "too 
close" to a smeared stellar residual. Thus, some of our detections fainter than my = 21 A 
undoubtedly arise from the confusing noise of bloomed stars (0BJ2 is far removed from 
any such noise). We use the population of noise objects detected in the reverse survey to 
set upper limits on the cumulative sky density of KBOs fainter than my = 27.4. Details of 
this calculation follow in the next section. 

The area searched in our deep survey is less than that of our shallow survey because 
of the shifting process. Areal losses ranged from 7-13% depending on the value of /i. To 
simplify the analysis, we adopt an average loss of -10%; the error introduced is negligible 
compared to Poisson uncertainties in the number of objects detected. Corrected for 
additional losses due to stellar /asteroidal residuals (-1%), our deep survey area equals 0.009 
square degrees = Ad. 



4.3. Cumulative Luminosity Function 

Figure || displays our estimates of the cumulative KBO sky density, S(< m^j), together 
with estimates made by various other groups. We emphasize that each survey's points 
represent estimates made independently of all other groups; i.e., survey areas have not been 
added.! Surveys conducted in V were included by assuming a solar color, V-R = 0.36, 
corresponding to a neutrally reflective KBO (red albedo pr = pv)- 

From our detection of OBJl, we independently estimate S(m/j < 25.14) = 98 ± 98(lo") 
objects/deg^. Combining this result with our detection of 0BJ2, which is the faintest KBO 
detected to date, we estimate T,{mR < 26.86) = 209 ± 149(lo") objects/deg^, where Poisson 



^The one exception to independence involves the points from Luu & Jewitt (1998), which presumably 
incorporate data from their previous surveys. 



-14- 



^|J, = -OJ A/^ = 0,0 A/x = +0.1 




Fig. 4. — Recombination images of 0BJ2. The central image is the best recombination 
image. Surrounding it are images recombined at adjacent points on the {fj,, 9) grid. Each 
panel to the right advances A/x — 0.1"/hr. Each panel towards the top of the page advances 
= 2°. 
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Fig. 5. — Comparison of the number of KBO candidates found in tiie deep forward survey 
and the number of noise objects detected in the reverse survey. At my < 27.4, the only object 
detected is 0BJ2, and no noise object is bright enough to confuse the identity of 0BJ2 as 
a true KBO. At my > 27.4, the number of KBO candidates never significantly exceeds the 
number of false alarms, and we can only compute upper limits on the sky density. 
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Fig. 6. — Independent estimates of the cumulative sky density of KBOs as made by various 
groups. Abbreviations for surveys are defined in the text and in the references. Upper hmits 
from this paper (CB99) are computed at the 99.99% ("4o"") confidence leveL Upper hmits 
from other surveys are pubhshed values at the 99% confidence level. 
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uncertainties have been added in quadrature. 

Upper limits are derived at fainter magnitudes as follows. We assume that in the 
forward survey, the occurrence of noise objects plus real KBOs is Poissonian. The expected 
mean number of forward survey candidates brighter than magnitude m in survey area 
Ad equals A^Noisc(< ^) + {v^d)^{< fn), where A''Noise(< ^) is the mean cumulative 
number of noise objects, and {i]Ad) is the efficiency-weighted survey area. We take 
N-^oise{< fn) = Nii{< m), where Nii{< m) is the cumulative number of reverse survey 
noise objects found. Given the number of forward survey candidates that we actually 
detected, Nf{< m), we ask what minimum value of S(< m) can be ruled out at the 99.99% 
confidence level (~"4(t" in Gaussian parlance): 



{Nn + (r/A,)S)^^ exp -{Nn + (r/A,)S) _ 

Npl ~ ^ ' ^ ^ 

an implicit equation for S where the magnitude dependence has been dropped for 
compactness. For reverse survey noise objects brighter than my (50%), we take 
{rjAd) = 0.7 X Ad- We do not calculate upper limits for my > my (50%), since our detection 
efficiency falls rapidly to zero past that magnitude (see Figure 1^). Upper limits on S 
computed using equation (H) are plotted in Figures ^ and |[ 

While Figure || summarizes the history of KBO surveys, quantitative results such as 
the slope of the luminosity function (or even the degree to which E resembles a single-slope 
power law) are better extracted from a fairer pooling of the data. To this end, we imagine 
the areas from all surveys as being combined into one giant frame over which the detection 
efficiency varies. At magnitude of a detected KBO, 



s(<"^.) = E , (3) 

k=l 

where rrij is the magnitude of the j^^ brightest KBO, rjk x Ak is the efficiency-weighted 
area of the k^^ survey, and n is the total number of surveys. Most surveys have published 
efficiency functions. Exceptions include the Mauna Kea-Cerro Tololo survey of Jewitt, Luu, 
& Chen (1996, hereafter JLC96), the Keck survey by Luu & Jewitt (1998, hereafter LJ98), 
the McGraw-Hill CCD survey by Luu & Jewitt (1988, hereafter LJ88), and the U.S. Naval 
Observatory survey by Levison & Duncan (1990, hereafter LD90). For data from JLC96, 
we assume i] behaves in a similar manner to that described in their companion paper I 
(Jewitt & Luu 1995, hereafter JL95); i.e., t] is assumed to fall linearly from 100% to 0% over 
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0.7 mag centered on published values of mji{50%). For data from LJ98, we obtained r] by 
consulting the lead author (Luu 1999). For the 2 remaining older surveys [which detected 
no KBOs, but which nonetheless contribute slightly to the total survey area in equation 
(^)], we adopted Heaviside step functions centered at uir = 24 (LJ88) and rriR = 22.14 
(LD90). We have verified that the conclusions of our paper are not affected by how we 
incorporate the latter 2 surveys. No photographic survey was included in the pool. 

Figure |^ displays the results of pooling datasets according to equation (H). To clarify 
the roles played by individual surveys, we pool an incrementally larger number of surveys 
in Figures ^ through |^c. In these and subsequent plots, magnitudes of individual points 
are identical to magnitudes rrii of individual KBOs. However, only points separated by 
at least ~0.2 mag are plotted; this represents a minor smoothing of the dataset, but is 
still preferable to imposing arbitrary bin boundaries. Error bars reflect Poisson counting 
statistics. 

Consider first Figure 0a, which incorporates data from Irwin, Tremaine, & Zytkow 
(1995), Jewitt, Luu, & Trujillo (1998), and Gladman et al. (1998) (hereafter 195, JLT98, 
and G98, respectively). These constitute the 3 surveys preferred by Gladman et al. (1998), 
excluding upper limit data. The points are well described by a power law, written in 
conventional notation as 

S(<m^) = 10«(^«-^o), (4) 

where slope a and reference magnitude mo are fitted parameters. A least-squares fit to 
these 3 surveys alone yields a = 0.70 ± 0.07, mo = 23.3 ± 0.1. These values coincide with 
those derived using a maximum likelihood analysis by Gladman et al. (1998); see their 
Figure 6c. We realize that least squares is not the preferred statistic for data whose errors 
are not Gaussian and which are correlated from point to point. However, the agreement 
between our result and G98's suggests that differences between slopes derived by various 
groups are due mainly to which surveys are kept and which are neglected, and not to the 
method of analysis. This will be borne out in what follows. 

In Figure 0b, we incorporate our survey (hereafter CB99) into the pool. The slope is 
lowered slightly to a = 0.66 ± 0.04, but the change is negligible over the range of observed 
magnitudes. On the basis of these 4 surveys alone, our data extend the a ~ 0.7 law to 
mn ^ 26.9. 

In Figure we fold in the Keck survey of LJ98. The observed faint end of the 
luminosity function is suppressed by the weight of their relatively sparsely populated fields. 
The luminosity function still resembles a single-slope power law, but the refitted slope is 
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Fig. 7. — Cumulative sky density obtained by pooling surveys according to equation (|^). 
From panels (a) through (c), successively more surveys are pooled, as indicated by the 
growing list of acronyms at the top of each panel. Data are fitted by least-squares, with 
fits from preceding panels plotted for comparison. In panel (a), the fitted slope a of the 
luminosity function is identical to that derived using a maximum likelihood analysis by 
Gladman et al. (1998). In panel (c), the slope a decreases significantly when data from 
LJ98 are included. 
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substantially shallower; a = 0.53 ± 0.02, coincident with the value given by LJ98. Though 
the LJ98 fields have a few times fewer objects at m/j ^ 26.5 compared to our CB99 fields, 
discrepancies are at the ~lcr level or less; uncertainties in our points (see Figure |^) are 
large because we detect only 2 objects. 

Finally, in Figure |], the remaining surveys by JL95, JLC96, LJ88, and LD90 are 
assimilated into the pool. The fitted luminosity function hardly changes; for this final pool, 
a = 0.52 ± 0.02 and tjiq = 23.5 ± 0.06. We note that shallow values for the fitted slope 
depend not only on surveys by JL95 and JLC96, but also on recent data from LJ98 (see 
Table 2). Some crude, mo del- dependent considerations of why values of a < 0.6 might be 
preferred are given in § [5.1| . 

The claimed Hubble Space Telescope (HST) detection of 2.5 x 10"^ objects/deg^ at 
rriR ~ 28.2 (Cochran et al. 1995) lies ~10 times above the steepest extrapolation, and 
~100 times above the extrapolation derived from all surveys combined. Brown, Kulkarni, 
& Liggett (1997) independently suggest on statistical grounds that the detections are 
erroneous. The increasing difficulty of reconciling the ground-based observations with the 
HST claim appears to support this suggestion. 

5. DISCUSSION 

5.1. Size, Surface Area, and Mass Distributions 

We consider a power-law differential size distribution with index q, dN{s) oc s~'^ds, 
where dN (s) is the number density of objects having diameters between s and s + ds. If 
all observed KBOs had the same albedo and were at the same heliocentric distance, the 
measured slope a of the cumulative luminosity function would imply a unique size index, 
viz. 

g = 5.02a + l. (5) 

This relation is straightforward to derive and is first given by 195. The assumption 
of uniform distance is not a bad one, since KBOs detected to date have present-day 
heliocentric distances between 30 and 50 AU; adopting a geometric mean distance for all 
objects mis-estimates sizes by at most a factor of 5/3. This is less than the possible factor 
of 4 uncertainty in size introduced by the unknown albedo, which might range from 4% 
(Comet Halley) to 60% (Pluto). JLT98 employ Monte Carlo models which incorporate 
more realistic distance distributions to extract the size index from the measured luminosity 
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Fig. 8. — Cumulative sky density obtained by pooling all surveys. The solid line is the 
power law fitted to all survey data. Other lines are fits from Figure ^, re-plotted here for 
comparison. Neither the HST datum nor Tombaugh's (1961) datum is included in any fit. 
Data from JL95 and JLC96 reinforce the shallow slope forced by data from LJ98. 
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function. Their best-fit q = 4.0 ±0.5 agrees with the value obtained by inserting their 
measured a = 0.58 ± 0.05 into equation (H). We shall use equation (||) to calculate q from 
a below, keeping in mind that such g's may be uncertain by ±0.5. 

Table § summarizes possible values of the sky density slope a and the size index q and 
their implications. Depending on which surveys are incorporated, q takes values from 3.6 
to 4.5. We compare these values to those of erosive disks in our Solar System. Main-belt 
asteroids are inferred to obey q ~ 3.3 in the diameter range 3-30 km (Durda, Greenberg, & 
Jedicke 1998). A value of g = 3.5 corresponds to a quasi-steady-state population for which 
catastrophic collisions grind as much mass per time into every size bin as they grind out, 
as first derived by Dohnanyi (1969). The derivation further assumes that critical specific 
energies for shattering and dispersal are independent of size.i This is a fair assumption for 
solid rocky targets smaller than ~10 km in diameter for which internal compression due to 
self-gravity is negligible. For asteroids greater than 30 km in diameter, there are significant 
deviations from the q ~ 3.5 law, with slopes ranging from g ^ 2 to 4.5 (Durda & Dermott 
1997). These deviations likely result from variations of the impact strength with size, as 
caused by self-gravitational effects (Durda, Greenberg, & Jedicke 1998; Melosh & Ryan 
1997). Saturn's ring particles crudely fit g ~ 3.3 in the size range of a few centimeters to 
a few meters based on Voyager radio occultation data (Marouf et al. 1983), though values 
between 2.8 and 4.0 cannot be completely ruled out (Cuzzi et al. 1984, Weidenschilling et 
al. 1984). Ring optical depths are sufficiently high that particles have suffered many erosive 
collisions over the age of the Solar System, so that their size distribution no longer purely 
reflects initial conditions (Borderies, Goldreich, & Tremaine 1984). 

For g > 3 (a > 0.4), surface areas (geometric optical depths) are dominated by the 
smallest bodies. All current estimates of g imply that this is the case for the Kuiper Belt. 

For g < 4 (a < 0.6), the total mass is dominated by the largest bodies. If we combine 
all surveys, we infer a size index g = 3.6 ±0.1. We use this g to make an order-of-magnitude 
estimate of the mass in the observable Kuiper Belt. Nearly all KBOs in the surveys we 
have considered have inferred diameters s 50 (0.04/p^)^/^ km.0 At limiting magnitude 
rrifi = 27 (the V-R adjusted magnitude above which false alarms prevent additional KBO 
detections in our deep survey), objects having s ^ 50 km can be seen out to distances of 48 
AU. For values of g and S(m/j < 27) derived by combining all survey data, the total mass 



^The critical specific energy for shattering, Q^, is defined as the energy per unit target mass required to 
produce a fragment with 50% the mass of the original target. It is smaller than Q^, the energy per unit 
target mass required to disperse such fragments to infinity (Melosh & Ryan 1997). 

i°The only exception is KBO K3, for which .s = 23km (Luu & Jewitt 1998). 
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Table 1. Properties of Detected KBOs 



Object 






e\deg) 


i?^(AU) 


z'^(deg) 


s'^(km) 






Method 


OBJl 


25.5 


3.83 


-1 


32.9 


4.5 


56 


100% 


0.0102 


Shallow 


0BJ2 


27.2 


2.92 





43.9 





46 


98% 


0.009 


Deep 



^Measured visual magnitude, uncertain by 0.3 (0.22) mag for OBJl (0BJ2). 

'^Measured proper motion amplitude and angle relative to ecliptic on CCD, respectively. For 
OBJl (0BJ2), uncertainties are 0.02 (0.05) "/hi and 1 (2.3) degrees. 

'^Inferred heliocentric distance and inclination, respectively, for an assumed circular orbit. 
^Inferred diameter, assuming a visual albedo of 4%. 
*^Detection efficiency and area searched, respectively. 



Table 2. Measured and Inferred 



a 



Source 



Implication 



0.52 All data (Fig. 8) 

0.53 Omit JL95, JLC96 (Fig. 7c) 

0.57 Omit LJ98 

0.66 Omit JL95, JLC96, LJ98 (Fig. 7b) 

0.70 Omit JL95, JLC96, LJ98, CB99 (Fig. 7a) 



3.6 Mass in Largest Bodies 

3.7 Mass in Largest Bodies 
3.9 Mass in Largest Bodies 
4.3 Mass in Smallest Bodies 
4.5 Mass in Smallest Bodies 



^Power-law slope of cumulative luminosity function; see equation (H). 

^Differential size index derived from q = 5.02 a+1, which assumes uniform albedo 
and distance. See 5|0| for discussion. 



of the Kuiper Belt out to 48 AU is 



M,e.(a< 48 AU) - 0.22 ^^^"^^ < ^ ^ "^^^ 



53deg-^ y VlO^degV V2gcm-3y 

Me . (6) 



0.04 \ / Sjyiax \ 



PR / V 2000 km. 



Here Akb is the sohd angle subtended by the Kuiper Belt (taken to extend ±15° in ecliptic 
latitude), p is the internal mass density of KBOs, and s^ax is the diameter of the largest 
body in the distribution (taken to be similar to Pluto). Our rough estimate of ~0.2 
is consistent with the upper limit of 1.3 M® within 50 AU derived by Hamid, Marsden, 
& Whipple (1968) using measured cometary orbits (see Weissman 1995 and MDL99 for a 
discussion of upper mass limits). Note that this model predicts the existence of ~10 more 
Pluto-sized objects in the nearby Kuiper Belt. 

The number of 1-10 km sized comet progenitors in the Kuiper Belt may be similarly 
estimated; 



. .TXN . in /S(mR < 27)\ / Akb \ /0.04\^ Vlkm\^-^ 
Nc,^.M < 48 AU) « 1.4 X 10'» (^j^ ) ) (— ) (— ) comets . 

(7) 

where Sc is the minimum diameter of a comet. Our order-of-magnitude estimate compares 
favorably with the population of ~7 x 10^ comets between 30 and 50 AU required to supply 
the rate of Jupiter-family comets (Levison & Duncan 1997) .0 

Omitting data from various surveys while preserving the same magnitude coverage in 
the luminosity function raises the inferred value of q and places most of the mass of the 
observable Kuiper Belt into the smallest objects. Gladman et al. (1998) do not incorporate 
data from JL95, JLC96, LJ98 and CB99. Their maximum likelihood analysis, which can 



^^Our calculation ignores the fact that some surveys observe ±90° away in ecliptic longitude from Neptune 
where Plutinos (KBOs in 3:2 resonance with Neptune) reach perihelion (Malhotra 1996). These surveys might 
be expected to find an unrepresentatively high sky density of KBOs. In fact these surveys (JL95, JLC96) 
find lower sky densities than other surveys; see § |4.3| and section 5.1 of G98. 

^^The scattered KBO disk has also been proposed as an alternative source of short-period comets. Duncan 
& Levison (1997) estimate that only 6 x 10* comets are required in the scattered disk to supply the observed 
rate. 
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and does assimilate upper limit data from Luu & Jewitt (1988) and Levison & Duncan 
(1990), concludes that the sky density slope a = 0.76. Inserting this value into equation (|) 
yields a size index q = 4.8. As a separate example of a shallow slope based on omission of 
data, a least-squares fit to the luminosity function which omits points from JL95, JLC96 
and LJ98, and which does not incorporate upper limit data, yields g = 4.3 (see Figure 
^). Both size indices would place most of the mass of the observable Kuiper Belt into the 
smallest objects. Since the size of the smallest object in the distribution is unconstrained, 
we cannot estimate the mass of the Kuiper Belt using these g's. However, for any g > 4, 
there always exists an Smin below which upper limits for the total cometary Belt mass 
within 50 AU (~1.3 M®) are violated. For values of g = 4.3 and T,{mR < 27) = 200deg~^ 
derived from Figure 0b, this minimum value for s^m is as large as 2 km, and only increases 
with increasing q. Explaining the existence of such lower cut-off sizes would be problematic. 

Our preferred size index, q = 3.6±0.1, is that of a Dohnanyi-like size distribution 
for objects having diameters between 50 and 500 km within 50 AU. However, whether 
the shape of this distribution results from a catastrophic collisional cascade as envisioned 
in Dohnanyi's (1969) scenario is questionable. The answer depends on impact strengths, 
relative velocities, and initial populations of KBOs, all of which are poorly constrained. 
For solid rocky bodies 50-500 km in diameter, critical specific energies for shattering and 
disruption are expected to increase strongly with size due to self-gravitational compression 
(Melosh & Ryan 1997, and references therein). The role of self-gravity is magnified yet 
further if bodies consist predominantly of weaker ices. Whatever their composition, we 
would not expect Dohnanyi's (1969) derivation to apply to objects as large as those 
observed, since the derivation assumes that impact strengths are independent of size. 
Relative velocities required for fragmentation and dispersal of solid rocky bodies ~100 km 
in size demand KBO eccentricities and inclinations exceeding 0.3; the actual history of 
the velocity dispersion is unknown. If KBOs consist of solid rock and relative velocities 
are sufficiently high for disruption and dispersal upon impact, we estimate that lifetimes 
against catastrophic dispersal of targets ~100 km in diameter exceed the age of the Solar 
System by a factor of ~150 if projectiles are drawn from the present-day Kuiper Belt. This 
estimate agrees with that of Stern (1995); see his Figure 2. Shaping the population of 
objects having sizes 50-500 km by catastrophic collisions would require a primordial Belt 
orders of magnitude more populous than what is observed today. 
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5.2. A Kuiper Cliff at 50 AU? 



To date not one member of the classical Kuiper Belt has been discovered beyond 50 
AU, despite observational advances in limiting magnitude and theoretical assurances that 
the region is dynamically stable. Gladman et al. (1998) have addressed this issue and 
concluded that the present sample of ~100 KBOs is marginally large enough to expect 
detection of such bodies. Here we confirm and elaborate upon their results. 

We assume the Kuiper Belt begins at an inner edge amin, and that the number density 
of objects (number per volume) decreases with distance a as a power law with index 



For a surface density (number per disk face area) appropriate to the minimum-mass outer 
solar nebula, the index /3 may plausibly take values of ~2-3, depending on how quickly 
random eccentricities e and inclinations i decay with heliocentric distance. In a field of 
limiting magnitude m, the sky density of objects (number per projected sky solid angle) 
located beyond distance a* is proportional to 



where 5^(0*) is the size of the smallest object which can just be seen at a*, and 
ci-,n{s) = a^^\J s/ Smia*) is the maximum distance out to which an object of size s can be seen. 
One immediate consequence of a Belt having distance and size indices considered here is 
that the faint end of the luminosity function is dominated by small nearby objects rather 
than large distant ones. Extending the limiting magnitude of a visual survey inherently 
achieves greater dynamic range in observable sizes than in distances because reflected 
fluxes decrease as s^/a^. The greater sensitivity to size is compounded by the shapes of 
the distributions; E oc s^~^a^~^ s~'^-^ aP'^ . An outer edge to such a Belt at 50 AU (a 
"Kuiper Chff " ) would not significantly flatten the slope of the luminosity function at faint 
magnitudes, a point which we shall justify more formally below. 

The fraction of objects located beyond a* is 



dN{s,a) oc a'/^s'^ds. 



(8) 





/ = 



E (a > ajnin) 



E (a > a*) 




1 + 




(10) 



where ^ — 2q + P — 5 (cf . G98) . The order-of-magnitude correction term is valid for 



-27- 



/3 < 3 and is small for surveys and distributions considered here. El The fraction / is thus 
insensitive to the limiting magnitude of the field. This insensitivity justifies our assertion 
that a Kuiper Cliff would not break the luminosity function at any particular magnitude. It 
also allows us to easily estimate how many detections beyond 50 AU we might expect. For 
flmm = 30 AU, a* = 50 AU, Smax = 2000 km, q = 3.6, and (3 = 3 (constant dispersion in e 
and i), the fraction of objects outside 50 AU is / ~ 8%. Decreasing the distance index /3 to 
2 increases / to 13%. In the present total sample of ~100 KBOs, we might therefore expect 
~10 to reside beyond 50 AU. Eight of these ten would be located between 50 and 70 AU. 

While these rough considerations do not convincingly implicate a Kuiper Cliff, they do 
argue more strongly against a sudden rise by factors of 3 or more in the surface number 
density between 50 and 70 AU (a nearby "Kuiper Wall"). Keeping the size distribution 
fixed and multiplying the surface density by 3 beyond 50 AU would demand that ~25% of 
all detected classical KBOs reside in such a wall, in contrast to the 0% found to date. 

Nonetheless, there are a number of ways the present lack of detections may still accord 
with a massive outer classical Belt. The size distribution of objects may change dramatically 
past 50 AU. For instance, it might be that only a few large objects exist between 50 and 
70 AU. This may plausibly be the result of runaway accretion unimpeded by the presence 
of Neptune. Alternatively, if only objects smaller than ~30 km populate the outer Belt, 
detecting them is a task better suited to occultation surveys than to searches relying on 
reflected light. 

6. SUMMARY 

Our main results follows. 

1. We discovered 2 new Kuiper Belt Objects in a single Keck LRIS field. One object 
at my = 25.5 was found by blinking individual frames. It lies at a heliocentric 
distance of i? ~ 33 AU and has a diameter s = 56 (0.04/py)^/^ km. Another object at 
my = 27.2 was discovered by blinking shifted and co-added frames. For this second 
object, ^ 44 AU and s = 46 [OM/pvY^'^ km. 

2. We pooled all surveys to construct the cumulative luminosity function from m^j = 20 
to 27 (Figure |). At the faintest observed magnitude, T,{mR ^ 26.9) = 53 ± 20(lo") 



i3por 7 « 5 and Smax — 2000 km, the correction term is less than 0.1 for surveys having Umiting magnitudes 
mn > 22. AU surveys used to construct our luminosity function satisfy the latter requirement. 
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objects/square degree. The best-fit slope is a = 0.52 ± 0.02, where 
logiQE(< mil) = a{mR — mo). Differences in a reported in the literatiirc are 
due mainly to which survey data are incorporated. Values of a > 0.6 require the 
omission of surveys by JL95, JLC96, and LJ98. 

3. Our KBO luminosity function is consistent with a power- law size distribution with 
differential size index q — 3.6 ±0.1 for objects having diameters 50-500 km within 
50 AU. The distribution is such that the smallest objects possess most of the surface 
area, but the largest bodies contain the bulk of the mass. By extrapolating outside the 
observed range of sizes, we estimate to order-of-magnitudc that 0.2M@ and 1 x 10^"^ 
comet progenitors lie between 30 and 50 AU. Though our estimated size index is 
that of a Dohnanyi-like distribution, the interpretation that catastrophic collisions 
are responsible is questionable. Impact strengths against catastrophic disruption 
and dispersal probably increase strongly with size for objects greater than ~10 km 
in diameter, whereas the derivation by Dohnanyi (1969) assumes impact strength 
to be independent of size. Lifetimes against catastrophic dispersal of KBOs having 
diameters 50-500 km exceed the age of the Solar System by at least 2 orders of 
magnitude in the present-day Belt, assuming bodies consist of competent rock. 

4. A greater than threefold rise in the surface density of the Kuiper Belt just beyond 50 
AU would imply that more than 25% of detected objects lie outside that distance, 
assuming objects are distributed similarly in size at all distances. The absence of 
detections past 50 AU in the present sample of ~100 KBOs argues against this 
picture. A massive outer Belt may still be possible if only a few large objects exist 
between 50 and 70 AU, or if only objects smaller than ~30 km exist in the outer Belt. 
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